arXiv:1507.00106v3 [quant-ph] 6Jul2015 


Event based simulation of an EPR-B experiment by 
local hidden variables: epr-simple and epr-clocked 

Richard D. Gill 
5 July, 2015 


Abstract 


In this note, I analyse the data generated by M. Fodje’s simulation programs epr-simple and epr-clocked 
using appropriate modified Bell-CHSH type inequalities: the Larsson detection loophole adjusted CHSH, and 
the Larsson-Gill coincidence loophole adjusted CHSH. The experimental efficiencies turn out to be ry « 81% 
and 7 « 55% respectively, and the observed value of CHSH is (of course) well within the adjusted bounds. 

Author address: R.D. Cill, Mathematical Institute, Leiden University, http://www.math.leidenuniv.nl/~gill 


Introduction 


Michel Fodje has written two event-based simulation programs of EPR-B experiments called “epr-simple” 
and “epr-clocked”. The programs are written in the Python programming language and are freely available at 
https://github.com/minkwe/epr-simple and https://github.com/minkwe/epr-clocked. Descriptions are given 
at https://github.com/minkwe/epr-simple/blob/master/README.md and https://github.com/minkwe/ 
epr-clocked/blob / master/RE ADME.md. 

The two programs use respectively the detection loophole (Pearle, 1970) and the coincidence loophole 
(Pascazio, 1986) to reproduce (to a good approximation) the singlet state correlation function and thereby to 
violate the CHSH inequality. See the appendix for further information on the models. 

In this note I study the experimental efficiency of the two models in the CHSH setting. Eor the detection 
loophole, the efficiency y is defined to be the minimum over all setting pairs and over all permutations of the 
two parties Alice and Bob, of the probability that Party 1 detects a particle given that Party 2 has detected 
a particle. Eor the coincidence loophole it is defined in a similar way: the efficiency 7 is defined to be the 
minimum over all setting pairs and over all permutations of the two parties Alice and Bob, of the probability 
that Party 1 has a detection which is paired to Party 2’s detection, given that Party 2 has detected a particle. 

If either loophole is present in the experiment, then the CHSH inequality is not applicable, or to be more 
precise, the statement that local hidden variables cannot violate CHSH is not true. I refer the reader to 
the survey paper Larsson (2014); arXiv eprint http://arxiv.org/abs/1407.0363. One needs to make further 
(untestable) assumptions such as the “fair sampling” hypothesis in order to deduce impossibility of local 
hidden variables from violation of CHSH. However, it is not difficult to modify CHSH to take account of the 
possibly differential “post-selection” of particle pairs which is allowed by these two loopholes. The result is 
two bounds, replacing the usual bound “2”: 4 / 77—2 for the detection loophole, and 6/7 — 4 for the coincidence 
loophole; see Larsson (2014) formula (38) and formulas (50), (51), (52). Note that when 77 = 1 the detection 
loophole bound equals the usual CHSH bound 2, but as 77 decreases from 1 , the bound increases above 2, 
at some point passing the best quantum mechanical prediction 2^/2 (the Tsirelson bound) and later even 
achieving the absolute bound 4. The bound is sharp: one can come up with local hidden variable models 
which exactly achieve the bound at the given efficiency. In particular, with 77 = 2/3 the detection loophole 
bound is 4, saying that it is possible for three of the four CHSH correlations to achieve their natural upper 
limit -|-1 and one of them its lower natural limit -1. 

The coincidence loophole bound is also attainable, and for the same value of the efficiency, worse. In particular, 
already with 7 = 3/4 one can attain three perfect correlations and one perfect anti-correlation. 
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epr-simple 


The programme epr-simple uses the detection loophole (Pearle, 1970) so as to simulate violation of the CHSH 
inequality in a local-realistic way. The simulated experiment can be characterised as a pulsed experiment. At 
each of a long sequence of discrete time moments, two new particles are created at the source, and dispatched 
to two detectors. At the detectors, time and time again, a new pair of random settings is generated. The 
two particles are measured according to the settings and the outcome is either -|-1, -1, or 0; the latter 
corresponding to “no detection”. 

If either particle is not detected, the pair is rejected. 

epr-simple only outputs some summary statistics for the accepted pairs. I added a few lines of code to 
the program so that it also outputs the “missing data”. Not being an expert in Python programming, my 
additional code is pretty simple. 

First of all, I reduced the total number of iterations to 10 million. The original code has 50 million, and 
this leads to memory problems on the (virtual) Linux Mint system which I use for Python work. Secondly, I 
added a code line “numpy.random.seed(1234)” in the block of code called “class Simulation(object)”, so that 
identical results are obtained every time I run the code. This means that others should be able to reproduce 
the numerical results which I analyse here, exactly. 

Finally, in the part of the code which outputs the simulation results for a test of the CHSH inequality, I 
added some lines to preserve the measurement outcomes in the case either measurement results in “zero” and 
then to cross-tabulate the results. 

By the way, for his test of CHSH, Michel Fodje (thinking of the polarization measurements in quantum 
optics) took the angles 0 and 45 degrees for Alice’s settings, and 22.5 and 67.5 for Bob. I have changed these 
to 0 and 90 for Alice, and 45 and 135 for Bob, as is appropriate for a spin-half experiment. 

for k,(i,j) in enumerateC [(a,b),(a,bp), (ap,b), (ap, bp)]); 

selO = (adeg==i) & (bdeg==j) # New variable 

sel = (adeg==i) & (bdeg==j) & (alice[;,l] != 0.0) & (bob[;,l] != 0.0) 

Ai = alicefsel, 1 ] 

AiO = alicefselO, 1 ] # New variable 

Bj = bob[sel, 1 ] 

BjO = bobfselO, 1 ] # New variable 

print " 7 .s; E( 7.5 . If , 7.5 . If ) , AB=7.+0.2f, qM=7.+0.2f" 7 . (DESIG[k],i, j, (Ai*Bj ) .meanO , 
-numpy.cos(numpy.radians (j-i))) 


npp = ((AiO 

== 

1 ) & 

(BjO == 

1 )) .sumO 



# 

New 

variable 

npO = ((AiO 

== 

1 ) & 

(BjO == 

0 ) ) . sumO 



# 

New 

variable 

npm = ((AiO 

== 

1 ) & 

(BjO == 

-1) ) . sumO 



# 

New 

variable 

nOp = ((AiO 

== 

0 ) & 

(BjO == 

1 )) .sumO 



# 

New 

variable 

nOO = ((AiO 

== 

0 ) & 

(BjO == 

0 ) ) . sumO 



# 

New 

variable 

nOm = ((AiO 

== 

0 ) & 

(BjO == 

-1) ) . sumO 



# 

New 

variable 

nmp = ((AiO 

== 

-1) & 

(BjO == 

= l)).sum() 



# 

New 

variable 

nmO = ((AiO 

== 

-1) & 

(BjO == 

= 0) ) . sumO 



# 

New 

variable 

nmm = ((AiO 

== 

-1) & 

(BjO == 

= -1)) .sumO 



# 

New 

variable 

print npp, i 

tipO, 

npm 


# Print 

out 

extra 

dcL i Oj 


print nOp, i 

nOO, 

nOm 


# Print 

out 

extra 

da i a 


print nmp, i 

nmO, 

nmm 


# Print 

out 

extra 

da i a 



CHSH.appendC (Ai*Bj ) .meanO ) 

QM.appendC -numpy.cos(nuiiipy.radians(j-i)) ) 


I ran epr-simple, redirecting the output to a text file called “data.txt”. In a text editor, I deleted all but 12 
lines of that file; the lines containing the numbers which are read into R, and then printed out by the R code 
below. 
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data <- as.matrix(read.table("data.txt")) 
colnames (data) <- NULL 
data [1:3, ] 


## 


[,1] 

[,2] 

[,3] 

## 

[1,] 

211 

330 

1284 

## 

[2,] 

351 

66 

325 

## 

[3,] 

1270 

339 

239 

data [4: 

6 , ] 




## 


[,1] 

[,2] 

[,3] 

## 

[1,] 

1304 

339 

201 

## 

[2,] 

315 

56 

336 

## 

[3,] 

250 

340 

1158 

data [7 : 

9, ] 




## 


[,1] 

[,2] 

[,3] 

## 

[1,] 

213 

324 

1313 

## 

[2,] 

319 

49 

348 

## 

[3,] 

1234 

311 

210 

data [10: 

:12, ] 




## 


[,1] 

[,2] 

[,3] 

## 

[1,] 

205 

305 

1290 

## 

[2,] 

314 

60 

319 

## 

[3,] 

1163 

341 

240 


dim(data) <- c(3, 2, 2, 3) 


Outcomes <- as.character(c (1 , 0, -1)) 

Settings <- as.character(c (1 , 2)) 

dims <- list(AliceOut = Outcomes, Alicein = Settings, Bobln = Settings, BobOut = Outcomes) 

dimnames (data) <- dims 

data <- aperm(data, c(l, 4, 2, 3)) 

data 


## 

, , Alicein = 

1, Bobln 

## 




## 

BobOut 



## 

AliceOut 1 

0 

-1 

## 

1 211 

330 

1284 

## 

0 351 

66 

325 

## 

-1 1270 

339 

239 

## 




## 

, , Alicein = 

2, Bobln 

## 




## 

BobOut 
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## 

AliceOut 

1 

0 

-1 

## 

1 

1304 

339 

201 

## 

0 

315 

56 

336 

## 

-1 

250 

340 

1158 

## 





## 

, , Alicein = 

1, Bobln 

## 





## 

BobOut 


## 

AliceOut 

1 

0 

-1 

## 

1 

213 

324 

1313 

## 

0 

319 

49 

348 

## 

-1 

1234 

311 

210 

## 





## 

, , Alicein = 

2, Bobln 

## 





## 

BobOut 


## 

AliceOut 

1 

0 

-1 

## 

1 

205 

305 

1290 

## 

0 

314 

60 

319 

## 

-1 

1163 

341 

240 


rho <- function(D) (D[l, 1] + D[3, 3] - D[l, 3] - D[3, 1])/(D[1, 1] + D[3, 3] + D[l, 3] + D[3, 1]) 
corrs <- matrixCO, 2, 2) 

ford in 1:2) {for (j in 1:2) corrs[i, j] <- rho(data[ , , i, j])} 
contrast <- c(-l, +1, -1, -1) 

S <- sumCcorrs * contrast) 
corrs 


## [, 1 ] [, 2 ] 
## [1,] -0.7003995 -0.7151515 
## [2,] 0.6903536 -0.6928916 

S ## observed value of CHSH 


## [1] 2.798796 


2 * sqrt(2) ## QM prediction (Tsirelson bound) 


## [1] 2.828427 

We see a nice violation of CHSH. However, a large number of particle pairs have been rejected. 

eta <- function(D) (D[l, 1] + D[l, 3] + D[3, 1] + D[3, 3])/sum(D[ , -2]) 
etap <- function(D) (D[l, 1] + D[l, 3] + D[3, 1] + D[3, 3])/sum(D[-2, ]) 

efficiency <- matrix(0, 2, 2) 

for(i in 1:2) {for (j in 1:2) efficiency[i, j] <- eta(data[ , , i, j])} 
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, i, j])> 


efficiencyp <- matrixCO, 2, 2) 

ford in 1:2) {for (j in 1:2) ef f iciencyp [i, j] <- etap(data[ , 


efficiency; efficiencyp 


## 


[,1] 

[,2] 

## 

[1,] 

0.8163043 

0.8166071 

## 

[2,] 

0.8173401 

0.8207307 

## 


[,1] 

[,2] 

## 

[1,] 

0.8178601 

0.8238558 

## 

[2,] 

0.8109688 

0.8177201 


etamin <- min (efficiency, efficiencyp) 
etamin 


## [1] 0.8109688 

We see that the minimum over all setting pairs, and over the two permutations of the set of two parties 
{Alice, Bob}, of the probability that Party 1 has an outcome given Party 2 has an outcome, is ry « 81%. 

The “correct” bound to the post-selected CHSH quantity S is not 2, but i/rj — 2 (Larsson, 2014). 

S # CHSH 

## [1] 2.798796 
4 / etamin - 2 # bound 

## [1] 2.932372 

We see that at this experimental efficiency, the local realism bound to CHSH is about 2.9, just above the 
Tsirelson bound 2\/2 ^ 2.8. The simulation generates results just below the bound: also at about 2.8. The 
observed value of S, the quantum mechanical prediction 2\/2 and the adjusted CHSH bound are all quite 
close together: the simulation model is pretty close to optimal. 


epr-clocked 

The program epr-clocked uses the coincidence loophole (Pascazio, 1986). Michel Fodje calls this a “clocked 
experiment” meaning that time is continuous, the times of detection of particles are random and unpredictable. 
(I would have preferred to reserve the word “clocked” as synonym for “pulsed”). Because Alice and Bob’s 
particles have different, random, delays (influenced by the detector settings which they meet), one cannot 
identify which particles were originally part of which particle pairs. Moreover, a small number of particles 
did not get detected at all, compounding this problem. 

The experimenter scans through the data looking for detections which are within some short time interval of 
one another. This is called the detection window. Unpaired detections are discarded. 

I ran the program, setting the spin to equal 0.5, and an experiment of duration 10 seconds. I let Alice and 
Bob use the settings for a CHSH experiment: Alice uses angles 0 and 90 degrees. Bob uses angles 45 and 
135 degrees. (As in epr-simple, Michel Fodje took the angles corresponding to a polarization experiment 
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instead of a spin experiment). I set the numpy random seed to the values “1234”, “2345”, and “3456” prior 
to running the source program, Alice’s station program, and Bob’s station programme respectively. This 
should make my results exactly reproducible ... but didn’t quite, because the 10 second duration of the 
experiment is ten seconds in “real time”. It therefore depends on queries by the program of the actual time in 
the real world outside the computer, and this process itself can take different lengths of time on each new 
run. However, the difference between the data obtained in different runs (with the same seed) should be 
negligeable: the total number of particle pairs will vary slightly, but their initial segments should coincide. 

In order to get a strictly reproducible simulation, I rewrote one section of the program “source.py”. Below is 
my replacement code. Instead of running for 10 seconds of real time, the code simply generates exactly 200 
000 emissions. 

def runCself, duration=10 . 0) : 

N = 200000 
n = 1 

print "Generating spin-7,0. If particle pairs" 7, (self.n/2.0) 
while n <= N: 
self.emit 0 
n = n + 1 

self.save ( 'SrcLeft.npy.gz' , numpy.array (self.left)) 
self.save (' SrcRight.npy.gz' , numpy.array (self.right)) 
print 

print "7od particles in ' SrcLeft. npy. gz ' " 7. (len(self . left) ) 
print "7od particles in ' SrcRight. npy. gz ' " 7o (len(self . right) ) 

The standard output gave me the following information: 

# No. of detected particles 

# Alice: 199994 

# Bob: 199993 

# Calculation of expectation values 

# Settings N_ah 

# 0 , 45 27416 

# 0, 135 27512 

# 90, 45 27345 

# 90, 135 27425 

# CHSH: <= 2.0, Sim: 2.790, QM: 2.828 

Notice the total number of detected particles on either side, and the total numbers of coincidences for each 
of the four setting pairs. The total number of coincidence pairs is a bit more than 100 thousand; the total 
number of detections on either side is almost 200 thousand. We have a rather poor experimental efficiency of 
about 55%. 

Npairs <- 27416 + 27512 + 27345 + 27425 
Nsingles <- 199994 
Npairs; Nsingles 


## [1] 109698 
## [1] 199994 
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gamma <- Npairs/Nsingles 
gamma 

## [1] 0.5485065 
6 / gamma - 4 


## [1] 6.938796 

The “correct” bound to the coincidence-selected CHSH quantity S is not 2. In fact, we do not know it exactly, 
but it is conjectured to be 6/7 — 4 (Larsson, 2014). Here, 7 is the effective efficiency of the experiment 
measured as the chance that a detected particle on one side of the experiment will be accepted as part of a 
coincidence pair. Notice that at 7 = 1 (full efficiency) the adjusted bound is equal to the usual bound 2, but 
that as it decreases from 100% the bound rapidly increases. At 7 = 3/4 it reaches its natural maximum of 4. 

At the observed efficiency of about 55%, the corrected CHSH bound in this experiment is close to 7, far 
above the natural and absolute bound 4. 

Because the efficiency is lower than in epr-simple, while the proper bound (adjusted CHSH) is higher, this 
experiment is a good deal worse than the previous one in terms of efficiency. It wouldn’t be too difficult, 
at this level of experimental efficiency, to tune parameters of this model so as to get the observed value of 
CHSH up to its natural maximum of 4. 

Incidentally, running epr-clocked many times, I experienced quite a few failures of the program “analyse.py” 
which is supposed to extract the coincidence pairs from the two data files. It seems that the Larsson 
algorithm for finding the pairs, which Fodje has adopted for this part of the data analysis, is failing in some 
circumstances. I could not find out what was the cause of this. 

Fortunately it is now rather easy to import Numpy (“numerical python”) binary data files into R using the 
package “RcppCNPy”. It should also not be difficult to find a suitable alternative to Larsson’s matching 
algorithm in the computer science literature and probably freely available in C-|— I- libraries. Algorithms 
written in C-|—I- can often easily be made available in R via “Repp”. Hence one could replace Michel Fodje’s 
“analyse.py” by one’s own data analysis script; this would also allow a “proper” computation of the efficiency 
7, taking the minimum over the efficiencies for each setting pair and both permutations of the two parties 
Alice and Bob. 


References 

R. D. Gill (2015), Pearle’s Hidden-Variable Model Revisited, http://arxiv.org/abs/1505.04431 

J.-A. Larsson (2014), Loopholes in Bell Inequality Tests of Local Realism, J. Phys. A 47 424003 (2014). 
http://arxiv.org/abs/1407.0363 

S. Pascazio (1986), Time and Bell-type inequalities, Phys. Letters A 118, 47-53. 

P.M. Pearle (1970), Hidden-Variable Example Based upon Data Rejection, Phys. Rev. D 2, 1418-1425. 


Appendix 

The appendices below contain mathematical formulas of the two simulation models. I have done my best to 
extract these faithfully from the original Python code and accompanying explanations by the author Michel 
Fodje. I have earlier published translations into the R language: http://rpubs.com/gillll09/epr-simple, 
http://rpubs.com/gillll09/epr-clocked-core, http://rpubs.com/gillll09/epr-clocked-full 
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Appendix: epr-simple 


Here is a little simulation experiment with (my interpretation in R of) epr-simple. I plot the simulated 
correlation function and also the acceptance rate. With an effective sample size of iV « 0.8 x 10®, the statistical 
error in simulated estimated correlation coefficiencts is roughly of size 0.001, well below the resolution of the 
graphs plotted below. Thus the small visible deviations from the theoretical negative cosine curve are for real. 
1 simplify the model by taking spin =1/2. Formulas are further simplified by a sign flip of all measurement 
outcomes, which by the symmetries of the model, does not change the observed data statistics. 

set. seed(9875) 

## For reproducibility. Replace integer seed by your own, or delete this line 
## and let your computer dream up one for you (system time + process ID). 

## Measurement angles for setting 'a'; directions in the equatorial plane 

angles <- seqCfrom = 0, to = 360, by = 1) * 2 * pi/360 

K <- length (angles) 

corrs <- numeric (K) ## Container for correlations 

Ns <- numeric (K) ## Container for number of states 

beta <-0*2* pi/360 ## Measurement direction 'b' fixed, in equatorial plane 

b <- c(cos(beta), sin(beta)) ## Measurement vector 'b' 

M <- 10~6 ## Size of "pre-ensemble" 

## Use the same, single sample of 'M' realizations of hidden 
## states for all measurement directions. This saves a lot of time, 

## and reduces variance when we look at *differences*. 

e <- runif(M, 0, 2*pi) 
ep <- e + pi 

U <- runif(M) 

p <- (sin(U * pi / 2)~2)/2 

## Loop through measurement vectors 'a' (except last = 360 degrees = first) 

for (i in 1 :(K - 1) ) { 
alpha <- angles[i] 
ca <- cos (alpha - e) 
cb <- cos (beta - ep) 

A <- ifelse(abs (ca) > p, sign(ca) , 0) 

B <- ifelse(abs (cb) > p, sign(cb) , 0) 

AB <- A*B 

good <- AB != 0 

corrs[i] <- mean (AB[good]) 

Ns[i] <- sum(good) 

} 


corrs[K] <- corrs [1] 
Ns[K] <- Ns[l] 


Correlation 


plotCangles * 180/pi, corrs, type = "1", col = "blue", main = "epr-simple" , 
xlab = "Angle (degrees)", ylab = "Correlation") 
linesCangles * 180/pi, - cos(angles), col = "black") 
legendCx = 0, y = 1.0, legend = c ("epr-simple" , "neg cosine"), 

text.col = cC'blue", "black"), Ity = 1, col = cC'blue", "black")) 

epr-simple 


o 


lO 

o 


o 

o 


o 

I 


o 

’T 

0 50 100 150 200 250 300 350 

Angle (degrees) 

plot(angles * 180/pi, Ns / M, type = "1", col = "blue", xlab = "Angle (degrees)", 
main = "Rate of detected particle pairs", ylim = c(0, 1)) 
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Rate of detected particle pairs 



The detection loophole model used here is very simple. There is a hidden variable E uniformly distributed in 
[027r]. Independently thereof, there is a second hidden variable P taking values in [0,1/2]. Its distribution is 
determined by the relation P = sin^((7r/2)[/)/2 where U is uniform on [0,1]. Alice and Bob’s measurement 
outcomes are signcos(i? —a) and sign cos(i? — /3 + 7r) respectively, if each of their particles is detected. Alice’s 
particle is detected if and only if abs(cos(i? — a)) > P and Bob’s if and only if abs(cos(i? — /3)) > P. 

Pearle (1970) characterized mathematically the set of all probability distributions of P which would give us 
the singlet correlations exactly (and for measurement directions in space, not just in the plane). He also picks 
out one particularly simple model in the class. His special choice has P = — 1 € [0,1] where V is 

uniform on [1,4], first expressed in this way by myself in 2014, see http://rpubs.com/gillll09/pearle and Gill 
(2015.) 

Below, the simulation is modified accordingly: just one line of code is altered. Now the experimental and 
theoretical curves are indistinguishable. 

# p <- (sindJ * pi / 2)~2)/2 ## epi—simple 

p <- 2/sqrt(l + 3 * U) - 1 ## Pearle (1970) model 


## Loop through measurement vectors 'a' (except last = 360 degrees = first) 


for (i in 1 :(K - 1) ) { 
alpha <- angles[i] 
ca <- cos (alpha - e) 
cb <- cos (beta - ep) 

A <- ifelse(abs (ca) > p, sign(ca) , 0) 

B <- ifelse(abs (cb) > p, sign(cb) , 0) 

AB <- A*B 

good <- AB != 0 

corrs[i] <- mean (AB[good]) 

Ns[i] <- sum(good) 
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Correlation 


corrs[K] <- corrs[l] 
Ns[K] <- Ns[l] 


plotCangles * 180/pi, corrs, type = "1", col = "blue", main = "Pearle", 
xlab = "Angle (degrees)", ylab = "Correlation") 
linesCangles * 180/pi, - cos(angles), col = "black") 
legendCx = 0, y = 1.0, legend = cC'Pearle", "neg cosine"), 

text.col = cC'blue", "black"), Ity = 1, col = cC'blue", "black")) 

Pearle 


p 


in 

o 


o 

o 


in 

o 

I 

o 

’T 

0 50 100 150 200 250 300 350 

Angle (degrees) 

plotCangles * 180/pi, Ns / M, type = "1", col = "blue", xlab = "Angle (degrees)", 
main = "Rate of detected particle pairs", ylim = c(0, 1)) 
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Rate of detected particle pairs 



Angle (degrees) 


Appendix: epr-clocked 

Here is the core part of (my interpretation in R of) epr-clocked: the emission of a pair of particles, supposing 
that both are detected and identified as belonging to a pair. However, if the detection times of the two 
particles are too far apart (larger than the so-called coincidence window), the pair is rejected. This would 
correspond to running the full epr-clocked model with very low emission rate. 

I have simplified the model by fixing spin =1/2. To save computer time and memory, I re-use the hidden 
variables E (“e”) and U from epr-simple. 

Again, I plot the simulated correlation function and also the acceptance rate. 

coincWindow <- 0.0004 

ts <- pi * 0.03 ## timescale 

asym <- 0.98 ## asymmetry parameter 

p <- 0.5 * sin(U * pi / 6)~2 

ml <- runif(M, asym, 1) ## small random jitter, left 

mr <- runif(M, asym, 1) ## small random jitter, right 

## Loop through measurement vectors 'a' (except last = 360 degrees = first) 

for (i in 1 :(K - 1) ) { 
alpha <- angles[i] 

Cl <- (0.5/pi) * cos(e - alpha) ## cos(e-a), left 
Cr <- (0.5/pi) * cos(ep - beta) ## cos(e-a), right 
tdl <- ts * pmax(ml * p - abs(Cl), 0) ## time delays, left 

tdr <- ts * pmax(mr * p - abs(Cr), 0) ## time delays, right 

A <- sign(Cl) ## measurement outcomes, left 
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B <- sign(Cr) ## measurement outcomes, right 

AB <- A * B ## product of outcomes 

good <- abs(tdl-tdr) < coincWindow 
corrs [i] <- mean (AB[good]) 

Ns[i] <- sum(good) 

} 

corrs[K] <- corrs [1] 

Ns[K] <- Ns[l] 

plotCangles * 180/pi, corrs, type = "1", col = "blue", main = "epr-clocked" , 
xlab = "Angle (degrees)", ylab = "Correlation") 
linesCangles * 180/pi, - cos(angles), col = "black") 
legendCx = 0, y = 1.0, legend = c ("epr-clocked" , "neg cosine"), 

text.col = cC'blue", "black"), Ity = 1, col = cC'blue", "black")) 

epr-clocked 


c 

g 

*-t—» 

o 

o 

O 



Angle (degrees) 


plot(angles * 180/pi, Ns / M, type = "1", col = "blue", xlab = "Angle (degrees)", 
main = "Rate of accepted particle pairs", ylim = c(0, 1)) 
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Rate of accepted particle pairs 



Angle (degrees) 

I now remove the “small random jitter”, effectively setting the “asymmetry parameter” to 1. I rescale time so 
that the model is finally described in terms of a couple of standard probability distributions and two arbitrary 
constants. 

coincWindow <- 0.034 

p <- 8 * p ## Now p lies in the interval [0, 1] 

## Loop through measurement vectors 'a' (except last = 360 degrees = first) 

for (i in 1 :(K - 1) ) { 
alpha <- angles[i] 

Cl <- cos(e - alpha) ## cos(e-a), left 

Cr <- - cos(e - beta) ## - cos(e-b), right 

tdl <- pmaxCp - 1.28 * abs(Cl), 0) ## time delays, left 

tdr <- pmaxCp - 1.28 * abs(Cr), 0) ## time delays, right 

A <- sign(Cl) ## measurement outcomes, left 

B <- sign(Cr) ## measurement outcomes, right 

AB <- A * B ## product of outcomes 

good <- abs(tdl-tdr) < coincWindow 

corrs[i] <- mean (AB[good]) 

Ns[i] <- sum(good) 

} 

corrs[K] <- corrs [1] 

Ns[K] <- Ns[l] 

plotCangles * 180/pi, corrs, type = "1", col = "blue", main = "epr-clocked" , 
xlab = "Angle (degrees)", ylab = "Correlation") 
lines(angles * 180/pi, - cos(angles), col = "black") 
legend(x = 0, y = 1.0, legend = c ("epr-clocked" , "neg cosine"). 
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"Rate of accepted particle pairs", ylim = c(0, 1)) 
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The results are almost identical. This allows us now to write down the core model in mathematics. The 
model is rather simple. There are just two hidden variables: a uniformly distributed angle E in [0, 27r] and 
independently thereof, a random number P = 4sin^(0) in [0,1] where 0 is uniformly distributed in [0,7r/6]. 

At Alice’s measurement device, where the setting is a, the measurement outcome is signcos(i? — a). At Bob’s 
measurement device, where the setting is P, the measurement outcome is signcos(A — (3 + tt). 

During measurement, the particles experience time delays. Alice’s particle’s time delay is max(P—1.28| cos{E— 
a) I, 0) and Bob’s is max(P — 1.28| cos{E — /3)|, 0). Notice that if a = /3 or if a = /3 + tt, the time delays of 
the two particles are identically equal to one another. 

Finally, the two detections are accepted as belonging to one particle pair if the difference between their two 
delay times is less than 0.034; i.e., if they are detected within the same time interval of length maximally 
0.034. 

The full “epr-clocked” model adds on top of this simple core, some further (relatively small) sources of 
noise, which do serve to smooth out the anomalous spike when the two particles are measured in the same 
direction. Moreover, in his simulations, Michel Fodje does not measure at fixed directions, but samples 
measurement directions uniformly at random in the circle. When he computes correlations, he has to bin 
measurement directions, resulting in another source of noise, again further smoothing out anomalous features 
in the observed correlation curve. 


Appendix: epr-clocked optimized 


Inspection of the data-sets generated by the last model turned up something very interesting: the coincidence 
window is so small, that for most measurement settings, a pair of detections is only accepted as a pair if both 
particles experience the same, zero, time delay! But this means that the two detections are accepted as a pair 
if and only if P < 1.28| cos{E — a) and P < 1.28| cos{E — /3)|, because only in this case are both time delays 
identically equal to zero. 

Now we are free to divide throughout in these inequalities by the constant 1.28 and absorb it into the random 
variable P, and we are free to pick a different distribution for P. So let us take P the same as in the Pearle 
model! Let us also reduce the size of the coincidence window to make it almost impossible for particles which 
experience non-zero time delays to become paired with their partners. 

coincWindow <- 0.000001 

p <- 2/sqrt(l + 3 * U) - 1 ## Pearle's choice 


## Loop through measurement vectors 'a' (except last = 360 degrees = first) 


for (i in 1 :(K - 1) ) { 
alpha <- angles[i] 

Cl <- cos(e - alpha) 

Cr <- - cos(e - beta) 
tdl <- pmaxCp - abs(Cl) 
tdr <- pmaxCp - abs(Cr) 
A <- sign(Cl) 

B <- sign(Cr) 

AB <- A * B 


## cos(e-a), left 
## - cos(e-b), right 
, 0) ## time delays, left 

, 0) ## time delays, right 

## measurement outcomes, 
## measurement outcomes, 
## product of outcomes 


good <- abs(tdl-tdr) < coincWindow 
corrs[i] <- mean (AB[good]) 

Ns [i] <- sum(good) 


left 

right 


corrs[K] <- corrs [1] 
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Correlation 


Ns[K] <- Ns[l] 


plotCangles * 180/pi, corrs, type = "1", col = "blue", 
main = "epr-clocked optimized", 
xlab = "Angle (degrees)", ylab = "Correlation") 
linesCangles * 180/pi, - cos(angles), col = "black") 
legend (x = 0, y = 1.0, legend = c ("epr-clock-opt" , "neg cosine"), 

text.col = cC'blue", "black"), Ity = 1, col = cC'blue", "black")) 

epr-clocked optimized 
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plotCangles * 180/pi, Ns / M, type = "1", col = "blue", xlab = "Angle (degrees)", 
main = "Rate of accepted particle pairs", ylim = c(0, 1)) 
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Conclusion: epr-clocked is actually an inferior version of the detection loophole model epr-simple, thinly 
disguised as a coincidence loophole model, but it is easy to vastly improve it. 
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